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ABSTRACT 

To examine a recently proposed hypothesis that silicon nanoparticlcs arc the source of extended red 
emission (ERE) in the interstellar medium, we performed a detailed modeling of the mean Galactic 
extinction in the presence of silicon nanoparticles. For this goal we used the appropriate optical constants 
of nanosized Si, essentially different from those of bulk Si due to quantum confinement. It was found that 
a dust mixture of silicon nanoparticles, bare graphite grains, silicate core-organic refractory mantle grains 
and three-layer silicate-water ice-organic refractory grains works well in explaining the extinction and, 
in addition, results in the acceptable fractions of UV/visible photons absorbed by silicon nanoparticles: 
0.071-0.081. Since these fractions barely agree with the fraction of UV/visible photons needed to excite 
the observed ERE, we conclude that the intrinsic photon conversion efficiency of the photoluminescence 
by silicon nanoparticles must be near 100%, if they are the source of the ERE. 

Subject headings: dust, extinction - ISM: general - methods: numerical 



1. INTRODUCTION 



& Berger 1996T) . The estimate of the amount of silicon 



Recently, extended red emission (ERE) has been de- 
tected in the diffuse galactic background and in high- 
galactic latitude cirrus clouds with surprisin gly high in- 
tensities (Gordon , Witt, & Friedmann 1998; pzomoru & 



Guhathakurta 1998). This emission, appearing in the 



500-800 nm spectral range as an unstructured broad fea- 
ture, is attributed to photoluminescence by some com- 
ponent of interstellar grains. The correlation between 
the ERE R-band intensity and H I-column de nsity of 
(1.43±0.31)1 0~ 29 ergs s" 1 A" ' 



sr 



1 H atom 1 ( |Gordon| 
et al. 1998 ) implies that, at a minimum, (10±3)% of 



all UV/visible photons in the diffuse interstellar radiation 
field are absorbed by the ERE producing grains, if each ab- 
sorption results in the emission of an ERE photon. This 
requires that the ERE-grains make a significant contribu- 
tion to the interstellar extinction in the UV/visible range, 
which further implies that they must consist of cosmically 
abundant elemental constituents in a form that permits 
extremely high efficiency photoluminescence. 

After a review of possible ERE-grain candidates, Witt, 
Gordon, & Furton (1998) and Ledoux et al. (1998) in- 
dependently proposed that silicon nanoparticles (SNPs) 
might be the source of ERE in the interstellar medium 
in as far as they produce photoluminescence spectra 
closely matching those observed in ERE sources and match 
the constraints regarding the required quantum efficiency 
and elemental abundance. It was discovered experimen- 
tally that the photoluminescence emission origin ates from 
nano c rystallinc Si region s in porous Si material flCanhair] 
1990 ; |Cole et al. 1992) an d from clusters of nanosized 
grains flTakagi et al. 1990t [Littau et al. 1994| ; |Tamii] 



needed in the form of SNPs, however, was based upon the 
average UV absorption cross section of crystalline silicon 
in bulk form. This must be considered as only a rough ap- 
proximation, because the remarkable optical properties of 
SNPs, clusters involving a few 10 2 silicon atoms each with 
a silicon suboxide (SiO^) surface covering, arise precisely 
because quantum confinement alters the band structure of 
the cluster material relative to that in the bulk material. 
Since quantum confinement is a function of cluster size, so 
are the resulting photoluminescence characteristics, and so 
are the dielectric functions of the cluster material. 

The current investigation was undertaken to answer sev- 
eral related questions: What are the UV/optical extinction 
characteristics of a separate component of SNPs? What 
fraction of the silicon believed to be depleted into solid 
grains is needed in the form of SNPs to explain the ob- 
served ERE intensities? What resulting characteristics of 
the interstellar extinction curve can be related to absorp- 
tion by SNPs and can thus be correlated with varying ERE 
intensities? To find answers on these questions, we made 
use of a powerful comput ational app roach based on the 
method of regularization ( [Zubko 1997 ), which was proven 
to be efficient in modeling of interstellar dust (Zubko, 
Krelowski, & Wegner 1996, 1998). 



2. OPTICAL PROPERTIES OF SI NANOPARTICLES 

The first step in our study was to choose the optical con- 
stants of SNPs. Witt et al. (1998) proposed that the in- 
terstellar SNPs should have diameters of 1.5-5 nm because 
only grains in this size range exhibit efficient photolumi- 
nescence. Evidently, the dielectric properties of such SNPs 
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should be different from those of a bulk material because 
of both the quantum size confinement effect and the in- 
creased role of surface states. These expect ations are fully 
confir med by avai lable experimental data ( Koshida et al] 
1993f| [rhcifi 1997| ). In particular, Theifi (1997) reported 
the dielectric function of nanosized silicon from porous 
samples with a range of porosities 60-80%, correspond- 
ing to the respective range of crystallite sizes of 5-2 nm. 
These results show that the optical properties of silicon 
nanostructures have a strong size dependence. However, 
an accurate derivation of the latter is not a trivial problem, 
because in each sili con sample t here is some size distribu- 
tion of crystallites ( [Theifi 1997j ). Note that Theifi' (1997) 
dielectric function covers the wavelengths 0.2-1 /im only. 
However, as we will see below, we especially need the data 
down to the Lyman limit (912 A). We therefore decided 
to use in our modeling the optical constants of nanosized 
silicon by Koshida et al. (1993) spanning the range 0.055- 
0.85 /im, which is sufficient for our problem. These con- 
stants correspond to the characteristic diameter of a SNP 
of around 3 nm. The optical constants of nanosized silicon 
by Koshida et al. (1993) and of bulk silicon taken from 
the compilation by Palik (1991) are displayed in Fig. [l| 
There is a considerable difference between the two sets 
which may also be seen in the extinction efficiencies of 
the respective spherical grains shown in Fig. 0. The first 
three cases correspond to the SNPs used in our modeling 
(see Section 3 for more details). Curves 3 and 4 demon- 
strate the so-called extrinsic size effect for small particles 
when the bulk optical constants are used. The peak at 
0.118 /jm is due to the SiC-2 shell and is present also in the 
spectra of bulk material, whereas the peaks at 0.153 /im 
and 0.138 /im (for 3 and 4, respectively) are surface plas- 
mons due to the silicon core, specific for small particles. 
The peak position depends on the volume fraction of Si. 
In contrast to the latter cases, curve 1 shows the intrin- 
sic size effect when the optical constants themselves vary 
with grain size. Note that Draine & Lee (1984) already 
computed the size-dependent dielectric function of small 
graphite grains. However, the size dependence was taken 
into account in their calculation only for the free-carrier 
part of the dielectric function, whereas the part due to the 
interband transitions was set to be the same as for bulk 
graphite. In contrast to Draine & Lee (1984), the optical 
constants of nanosized silicon used in the present study 
explicitly reflect the drastic changes in the band structure 
and interband transitions of the SNPs in comparison with 
those of bulk silicon. 



3. INTERSTELLAR EXTINCTION WITH SI 
NANOPARTICLES 

The dust models investigated for the present study were 
expected to satisfy three independent constraints simul- 
taneously: 1) fit the mean Galactic extinction curve as 
closely as possible over the entire near-IR to far-UV range; 
2) consume depleted heavy elements consistent with the 
recently proposed reference abundances for the interstel- 
lar medium (about 2/3 of solar values, see e.g., Snow & 
Witt 1996); and 3) require the ERE-producing SNPs to 
absorb at least 0.10±0.03 of the UV/visible photons in 
the interstellar radiation field, as suggested by the obser- 
vations of Gordon et al. (1998). The problem was solved 



by using the approach and computer program based on the 
method of regularization (Zubko 1997; Zubko et al. 1998) 
which allows one to find the size distributions of the dust 
components in the most general form. The size distribu- 
tions may be subject to the cosmic abundance and mass 
fraction constraints. Following Zubko et al. (1998), we 
treated the mixtures of spherical bare, multilayer and/or 
composite grains, whose extinction efficiencies were calcu- 
lated either with the standard Mie theory for bare and 
composite grains (Bohren & Huffman 1983) or with its ex- 
tension to multilayer particles. For composite grains the 
effective dielectric function was first calculated using the 
effective medium approach by Stognicnko, Hcnning, & Os- 
senkopf (1995). 

The fraction of UV/visible photons absorbed by SNPs, 
P, was calculated with the formula 



,SNP )r SN P/ ISRF dA 



/o°oS m (l-«A)r A ir F dA 



(1) 



where a\ (a| NP ) is the albedo of the entire dust mixture 
(SNPs), t\ (r| INP ) is the extinction of the entire dust mix- 
ture (SNPs), /i SRF is the intensity of the interstellar ra- 
diation field at the Earth in units of photons cm~ 2 s _1 
A -1 . The a\ was taken from Gordon et al. (1998), t\, the 
mean Galactic extinction (i?y=3.1) from Cardelli, Clay- 
ton, & Mathis (1989) and if RF from Mathis, Mezger, & 
Panagia (1983). The values of a^ NP and t| np were taken 
from the models. 

We based our calculations on the G, M and O types of 
models introduced by Zubko et al. (1998). A G model 
contains the constituents usual for Greenberg's models: 
the bare graphite grains and the silicate core-organic re- 
fractory grains (e.g. Li & Greenberg 1997). An M model 
is a mixture of the bare graphite and silicate grains to- 
gether with the porous composite grains, resembling recent 
Mathis' (1996) models. The constituents of an O model 
are the bare graphite and silicate grains together with the 
multilayer grains consisting of a silicate core coated by re- 
fractory organic and water-ice shells. 

Initially, we tried to model the extinction using the 
above models to which we added SNPs of 3.0 nm diame- 
ter covered by a Si02 mantle. We varied the thickness of 
the mantle and found that the best results were achieved 
when the mantle is quite thin. However, all the models 
resulted in quite low (5, as a rule not exceeding 0.03-0.04, 
and, as a rule, consumed the solar fractional amount of 
silicon, 35 ppm (atoms per 10 6 H atoms). For the next 
step we introduced spherical porous clusters of SNPs in- 
stead of isolated SNPs. This resulted in an increase of f3, 
typically by 0.003-0.005 (the highest values correspond to 
70% porosity). Then we tried combining the G, M, and 
O models and found that the only mixture that results 
in acceptable (3 (> 0.07) is a combination of G and O 
models or the GO model. It consists of four components: 
bare graphite grains, silicate core-organic refractory man- 
tle grains, three-layer silicate- water ice-organic refractory 
grains and, finally, the SNPs in the form of either sili- 
con core-Si02 mantle grains or 70%-porous aggregates of 
silicon grains containing a small amount of Si02- The 
GO model requires a mass fraction of 7% of silicon grains. 
It results in /3=0.071 and 0.081 for the models with iso- 
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lated and clustered SNPs, respectively. Note, however, 
that these models still consume much carbon, 190 ppm, 
and a solar fractional amount of silicon. The GO model 
seems more natural than G and O models separately, be- 
cause when the dust grains in the interstellar medium do 
evolve through the formation and processing of mantles, 
we should expect the contributions into extinction of var- 
ious intermediate stages. 

The main results of our calculations are presented in 
Fig. | and Tables [|-§ First of all, we see that the mod- 
els with bulk silicon exhibit rather strong features in the 
UV (see Fig. |J) which so far have not been observed. 
For this case we adopted the volume fractions: <5si=0.297 
and <Sgio 2 =0.703 corresponding to the initial assumption 
of Witt el al. (1998) about an equal amount of silicon and 
oxygen atoms in a SNPs. By contrast, the contribution of 
the SNPs to the extinction is featureless when the optical 
constants of nanosized silicon and a very thin Si02 mantle 
are used. As a result, such models fit the mean Galac- 
tic extinction much better. Practically all of the models 
require much carbon: 180-190 ppm, except O models for 
which C/H=100 ppm, which is about optimal. As found 
in previous models by Zubko et al. (1998), more oxygen 
is needed for the O models. The amount of silicon needed 
requires solar relative abundances. A remarkable feature 
of the GO models in contrast to the other models is the 
extremely small amount of magnesium and iron consumed, 
about 17 ppm, which is quite close to the respective opti- 
mum values by Snow & Witt (1996). In addition to pro- 
viding the highest (3 values, the GO models also show the 
closest ht to the observed mean Galactic extinction curve. 

A result found to be generally true for all of our models 
is that (3 is approximately equal to the mass fraction incor- 
porated in the SNPs (see Table ||) . This result is largely 
independent of whether the SNPs are included as a single- 
size component of a=1.5 nm, an a -3 distribution of porous 
SNP clusters in the range 1.5-50 nm, or, as pointed out by 
the referee, if SNPs are included as the small-size tail of an 
a -3.5 distribution, as long as their mass remains constant. 
Since the cosmic abundance of Si limits this mass fraction 
severely and approximately 50% of the available Si is al- 



ready incorporated in the SNPs in the GO models, one 
important conclusion to draw is that if ERE is produced 
by SNPs, the intrinsic photon conversion efficiency must 
be extremely high, essentially near 100%. 

One reason why the (i values from Table || are as small 
as they are is that the number of photons in the assumed 
radiation field is highest in the 0.35-0.55 /im wavelength 
range, where the absorption contribution of the SNPs is 
practically zero. One would obtain larger values of (3, if 
the interstellar radiation field contained increased num- 
bers of UV photons. If one can make a case for significant 
differences in the spectral energy distribution of the in- 
terstellar radiation field in different regions of the Galaxy, 
one should be able to predict corresponding variations in 
the relative brightness of the ERE and the visible diffuse 
galactic light (DGL). Such variations should be apparent 
in results from ongoing surveys which map the DGL and 
the ERE over a w ide range of galactic longitudes (e.g., 
Gordon et al. 1998| ). 

We conclude that in general the results of our models are 
consistent with the hypothesis that SNPs are the source of 
ERE in the interstellar medium. Their extinction contri- 
bution occurs primarily at wavelengths shorter than 300 
nm and is not expected to introduce discernable structure 
into the interstellar extinction curve, provided one takes 
into account the size dependence of the optical constants 
for nanostructured silicon. The intrinsic photon conver- 
sion efficiency of such SNPs must be near 100%, if cosmic 
abundance constraints are to be met. Spherical clusters 
of SNPs with cluster radii up to 50 nm provide a slightly 
better match to the ERE requirement than do single SNPs 
with the same total mass fraction. It would be of great in- 
terest to predict possible spectral features due to SNPs in 
the infrared, as well as the contribution of SNPs in the 
total thermal emission of interstellar dust. However, we 
cannot do this because of the lack of optical constants of 
SNPs in the infrared and because of lacking information 
about the size distribution of SNPs. 

We thank Prof. N. Koshida for providing us the elec- 
tronic table of the optical constants of nanosized silicon. 



REFERENCES 



Bohren, C.F., & Huffman, D.R. 1983, Absorption and Scattering of 

Light by Small Particles (New York: Wiley-Interscience) 
Canham, L.T. 1990, Applied Physics Letters, 57, 1046 
Cole, M.W., Harvey, J.F., Lux, R.A., et al. 1992, Applied Physics 
Letters 60, 2800 

Cardelli, J. A., Clayton, G.C., & Mathis, J.S. 1989, ApJ, 345, 245 
Draine, B.T., & Lee, H.M. 1984, ApJ, 285, 89 
Gordon, K.D., Calzetti, D., & Witt, A.N. 1997, ApJ, 487, 625 
Gordon, K.D., Witt, A.N., & Friedmann, B.C., 1998, ApJ, 498, 522 
Koshida, N., Koyama, H., Suda, Y., et al. 1993, Applied Physics 

Letters, 63, 2774 
Ledoux, G., et al. 1998, A&A, 333, L39 
Li, A., & Greenberg, J.M., A&A, 323, 566 

Littau, K.A., et al. 1994, Journal of Physical Chemistry, 97, 1224 



Mathis, J.S., ApJ, 472, 643 

Mathis, J.S., Mezger, P.G., & Panagia, N. 1983, A&A, 128, 212 
Palik, E.D. 1991, Handbook of Optical Constants of Solids (Boston: 

Academic Press) 
Snow, T.P., & Witt, A.N. 1996, ApJ, 468, L65 

Stognienko, R., Henning, Th., & Ossenkopf, V. 1995, A&A, 296, 797 

Szomoru, A., & Guhathakurta, P. 1998, ApJ, 494, L93 

Takagi, H., et al. 1990, Applied Physics Letters, 56, 2379 

Tamir, S., & Berger, S. 1996, Thin Solid Films, 276, 108 

Their], W. 1997, Surface Science Reports, 29, 91 

Witt, A.N., Gordon, K.D., & Furton, D.G. 1998, ApJ, 501, Llll 

Zubko, V.G, 1997, MNRAS, 289, 305 

Zubko, V.G, Krelowski, J., & Wegner, W. 1996, MNRAS, 283, 577 
Zubko, V.G, Krelowski, J., & Wegner, W. 1998, MNRAS, 294, 548 



4 



Silicon nanoparticles and interstellar extinction 



Table 1 

The parameters of the GO models. 



D a 




Components 


f c 

J mass 


C d 


Si d 


O d 


J.VJ.g, 


Fc d 


n 


i 




100.0 


190 


35 


180 


17 


17 






Graphite 


22.1 


143 




















/inn 
4U.U 


go 
OO 


Iz 


QQ 
OO 


Iz 


Iz 






Sil-Ice-ORR 


30.9 


14 


5 


97 


5 


5 






Si-SiG-2 


7.0 





18 











n 


c 




100.0 


190 


35 


180 


17 


17 






Graphite 


22.1 


143 


















Sil-ORR 


40.0 


33 


12 


83 


12 


12 






Sil-Ice-ORR 


30.9 


14 


5 


96 


5 


5 






Si-Si02-voids 


7.0 





18 


1 








b 


i 




100.0 


190 


29 


180 


17 


17 






Graphite 


23.5 


145 


















Sil-ORR 


40.0 


32 


12 


80 


12 


12 






Sil-Ice-ORR 


29.5 


13 


5 


88 


5 


5 






Si-SiQ 2 


7.0 





12 


12 









a The optical constants of nanosized (n) or bulk (b) silicon. 
b The type of a SNP: isolated (i) or clustered (c). 
c The mass fractions of the dust components (/ maS s)- 

d The amounts of chemical elements in dust in ppm: atoms per 10 6 H atoms. 



Table 2 

The parameters of the models. 



Model 


D 


T 


A* 


/SNP 




GO 


n 




0.284 


0.070 


0.071 


GO 


n 




0.289 


0.070 


0.081 


GO 


b 




0.406 


0.070 


0.054 


G 


n 




0.297 


0.030 


0.033 


G 


n 




0.273 


0.030 


0.038 


G 


b 




0.315 


0.030 


0.026 


M 


n 




0.352 


0.020 


0.021 


M 


n 




0.353 


0.020 


0.024 


M 


b 




0.357 


0.020 


0.018 


O 


n 




0.373 


0.033 


0.037 


O 


n 




0.418 


0.033 


0.042 





b 




0.427 


0.033 


0.029 



a The mass fractions (/snp) and the fractions of UV/visible photons absorbed by SNPs ((3). 

b The quality of the fit is characterized by the normalized measure of incompatibility fi (Zubko et al. 1996). 
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Fig. 1. — The optical constants of bulk and nanosized silicon 
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Fig. 2. — The extinction efficiency per grain radius, Qext/a, of the SNPs and graphite grains. The optical constants of both nanosized 
(NS) and bulk (BS) silicon were used. Shown are the cases: 1: NS (0.99) core-Si0 2 (0.01) mantle, 2: NS(0.29), Si0 2 (0.01), voids (0.7) in 
a composite grain, 3: BS (0.297) corc-SiC>2 (0.703) mantle, 4: BS (0.700) core-Si02 (0.300) mantle, and 5: graphite grains. The respective 
volume fractions arc noted in brackets. 



6 Silicon nanoparticles and interstellar extinction 




I 10 100 1000 11 10 987654321 




1 10 100 1000 11 10 9 87654321 

grain radius, nm 1/^, urn" 1 



Fig. 3. — Size distributions of dust grains (left panels) and respective extinctiotLcurves (right panels), fitting the Galactic mean extinction 
curve (Rv=3.1). Shown are the GO-type dust models. See the text and Tables fU-felfor more details. 



